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Abstract 

Recent observations of neutrino oscillations imply non-zero neutrino masses and 
flavor violation in the lepton sector, most economically explained by the seesaw 
mechanism. Within the context of supersymmetry, lepton flavor violation (LFV) 
among the neutrinos can be communicated by renormalization group flow to the 
sleptons and from there to the charged leptons. We show that LFV can appear 
in the couplings of the neutral Higgs bosons, an effect that is strongly enhanced 
at large tan (3. In particular, we calculate the branching fraction for r — ► 3/x and 
fi — ► 3e mediated by Higgs and find that they can be as large as 10 and 5 x 10 -14 
respectively. These mode, along with B° — » fiu, can provide important evidence for 
supersymmetry before direct discovery of supersymmetric partners occurs. Along 
with t — > fjL'y and fx — > &y, they can also provide key insights into the form of the 
neutrino Yukawa mass matrix. 



Over the last several years, evidence from a number of experiments, notably Su- 
per Kamiokande and SNO, has pointed conclusively to the existence of neutrino oscilla- 
tions in atmospheric Q and solar neutrinos |2| and, by implication, to non-zero neutrino 
masses. Within the context of the Standard Model (SM), the most attractive explanation 
for the observed neutrino masses is the "seesaw" mechanism 0: right-handed neutri- 
nos are introduced in order to couple with left-handed neutrinos through SU(2) xU(l)- 
violating Dirac mass terms, tud, while also receiving large, SU(2) xU(l)-invariant Ma- 
jorana masses, Mr. The resulting spectrum consists of heavy neutrinos with masses 
~ Mr which are primarily right-handed, and neutrinos with extremely small masses 
m v ~ /Mr which are primarily left-handed. The seesaw mechanism also has the 
virtue that it fits elegantly inside a grand unified theory (GUT), such as SO (10). 

Within SO (10), the Dirac neutrino masses are predicted to be of order the correspond- 
ing up-quark masses; for example, [ra D ) VT would be roughly 100 to 200 GeV. Atmospheric 
neutrino data favors a u T mass of about 0.04 eV Thus one finds a right-handed Ma- 
jorana mass, Mr, of order 10 14 GeV, several orders of magnitude below the Planck scale 
(M Pl ). 

Majorana neutrino masses and neutrino oscillations imply lepton flavor violation 
(LFV). But within the SM, flavor violation in charged lepton processes is necessarily 
generated by irrelevant operators and is therefore suppressed by powers of 1/Mr (or 
equivalently powers of m„). Within supersymmetric extensions of the SM, however, this 
is no longer true. In the minimal supersymmetric standard model (MSSM) (which we 
will henceforth take to be augmented with heavy right-handed neutrinos, Ur) LFV can 
be communicated directly from vr to the sleptons by relevant operators and from there 
to the charged leptons. LFV is then suppressed by powers of 1/M SU sy instead of 1/Mr, 
with Msusy *C Mr. The initial communication is done most economically through 
renormalization group flow of the slepton mass matrices at energies between Mpi and 
Mr. Though the scale Mr is far above the weak scale, the presence of the Vr at scales 
above Mr leaves a (non-decoupling) imprint on the mass matrices of the sleptons which 
is preserved down to the weak scale. This effect has been used to predict large branching 
fractions for r — ► /ry and /i — > cy within the MSSM [0, 0, 0. 

In this letter we demonstrate a new way in which the imprint of LFV on the slepton 
mass matrices can be communicated to charged leptons through the exchange of Higgs 
bosons, providing the possibility of new and observable flavor violation in the leptonic 
sector. We will demonstrate that the decay r — > 3/i is a particularly sensitive probe of 
LFV at large tan (3 (the ratio (H u ) / (Ha)), with a branching fraction scaling as tan 6 (3. The 
process — > 3e can also proceed via Higgs exchange with the same tan 6 /3-dependence; 
though the branching fraction is suppressed by the small electron Yukawa coupling, it is 
still within the grasp of experiment. We will concentrate our discussion on the r decay 
and come back to the \i decay in the final discussion. 

At the end we will also briefly discuss the related processes r — > fifie and r — > /iee. 
We find no new enhancement for the decay r — > 3e. These decay modes are (like r — > fj,j) 
only logarithmically sensitive to Mr, but (unlike r — > fj,j) do not decouple in the large 
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mi mass limit. Rather, these modes decouple in the limit that the pseudoscalar Higgs 
boson becomes heavy, uia — > oo, thus providing complementary information on the 
super symmetric (SUSY) spectrum. 

Flavor Violation among the Sleptons. In the leptonic sector, we begin with a 
Lagrangian: 

- C = E R Y E L L H d + V R Y V L L + \v T R M R v R (1) 

where E R , Ll and v R represent 3x1 matrices in flavor space of right-handed charged 
leptons, left-handed lepton doublets and right-handed neutrinos, and Ye, Y v and M R are 
3x3 matrices in flavor space; for example, E R = (e R , fi R ,r R ) T . This Lagrangian clearly 
violates both family and total lepton number due to the presence of the v R Majorana 
mass term. We can choose to work in a basis in which both Ye and M R have been 
diagonalized, but Y v remains an arbitrary, complex matrix. 

Within the SM, 0(1) flavor violation in the neutrinos does not translate into appre- 
ciable flavor violation in the charged lepton sector due to 1/Mr suppressions. But this 
is not true in the slepton sector. The SUSY-breaking slepton masses are unprotected 
by chiral symmetries and are therefore sensitive to physics at all mass scales between 
mi and the scale, M, at which SUSY-breaking is communicated to the visible sector, 
assuming M > M R . This can be seen by examining the renormalization group equation 
for m| at scales above M R : 



+ 



/ MSSM 
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m\Y}Y v + YlY v m\ + 2(Y}m* H Y v + m 2 H Y} Y >> + A t A » 



where the first term represents the (L-conserving) terms present in the usual MSSM at 
scales below Mr. Because Y v is off-diagonal, it will generate flavor-mixing in the slepton 
mass matrix. We can solve this equation approximately for the flavor-mixing piece: 

- {^urX), + 2 (ma\) ( 3) 

where mo is a common scalar mass evaluated at the scale Q = M, and % ^ j. If we 
further assume that the A-terms are proportional to Yukawa matrices, then: 

( Am lL - £ { Y ' Y ") H (4) 



£ = T^Zl (6 + 2a)m . (5) 



where 

\og(M/M R ) 
16tt 2 

and a is 0(1). In the simplest SUSY-breaking scenarios, gravity plays the role of mes- 
senger and M = Mpi, so that the logarithm in Eq. (H) is roughly 10. 
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What does experiment tell us about the values of these matrices? Global fits to 
neutrino data favor large mixing between the and u T , and also between v e and |4[]. 
We will consider the following form for m v which provides an excellent fit to existing 
neutrino data and can be motivated by theory 




m v oc e 1 1 (6) 



where e is a small parameter ~ 0.1. If we further assume that Mr is an identity matrix, 
then YjY v will also have the form of Eq. Another interesting possibility is provided 
by GUT models with lopsided mass matrices for charged leptons 0; such models have 
{Ye)^2 — (Ye)33 and lead to a light neutrino mass matrix as in Eq. (H) with (Y u ) 3 2 — 
(1^)33 ~ y t , where y t is the top Yukawa coupling. In either case, the YjY v has 0(1) flavor 
violation in the v T —v )li sector if Mr ~ 10 14 GeV. 

Higgs-Mediated Flavor Violation. Unlike the SM, the MSSM is not protected 
against the possibility of flavor- changing neutral currents (FCNCs) mediated by neutral 
Higgs bosons. Though the MSSM is a type-II two-Higgs doublet model at tree level, 
this structure is not protected by any symmetry. In particular, the presence of a non- 
zero //-term, coupled with SUSY-breaking, is enough to induce non-holomorphic Yukawa 
interactions for the quarks and leptons. In the quark sector this was discovered by Hall, 
Rattazzi and Sarid | T0| ; terms of the form Q l urH\ and Q L d R Hl were found, the latter 



providing a significant correction to the 6-quark mass at large tan (3. We have shown in a 



previous letter [ITJ that these terms allow the neutral Higgs bosons to mediate FCNCs, 
in particular B — > ///i. There we argued that branching fractions predicted at large 
tan p can be probed at Run II of the Tevatron. (For recent analyses of Higgs-mediated 
B^ufi, see Ref. O.) 



The two leading diagrams considered in Refs. |1(J, [UJ as a source for non-holomorphic 
quark couplings are not present in the leptonic sector since they involve gluinos and top 
squarks inside the loops. However there are additional diagrams which are present in 
the leptonic sector JTJ] involving loops of sleptons and charginos/neutralinos; a subset of 
these is shown in Fig. 1. Given a source of non-holomorphic couplings and LFV among 
the sleptons, Higgs-mediated LFV is unavoidable. 



For those familiar with the quark sector FCNC calculation of Ref. |L1] , we note that 
the present calculation is not very different. In particular, the contributions considered 
here will be similar in most ways to those in Ref. |11] which were generated by a squark 



mixing insertion in the gluino-sbottom loop diagram. 

We begin by writing the effective Lagrangian for the couplings of the charged leptons 
to the neutral Higgs fields: 

- C = E R Y E E L H° d + E R Y E (e x l + e 2 YjY v ) E L H° U * + h.c. (7) 

The first term is the usual Yukawa coupling, while the second term arises from the 
non-holomorphic loop corrections. LFV results from our inability to simultaneously 
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Figure 1: Diagrams that contribute to e 2 . The crosses on the internal slepton lines represent 
LFV mass insertions due to loops of vr. 

diagonalize Ye and the 6 2 YeYJY u termQ, just as in Ref. [pi]] ; as e 2 — > 0, LFV in the Higgs 
sector will disappear. 

The diagrams which contribute to e 2 are shown in Fig. 1. Each diagram contains a 
single insertion of Am| which introduces LFV into the process. Without this insertion, 
these diagrams would have a trivial flavor structure and would not contribute to e 2 or to 
LFV. But the Am| insertion introduces a YJY V into the diagram, yielding a contribution 
to e 2 . We can approximate the contributions of the diagrams in Fig. 1 by inserting a single 
Am- mass insertion onto each of the internal El lines. We will also treat the higgsinos 
and gauginos as approximate mass eigenstates. For diagram 1(a), the contribution to e 2 



is 

a' 



e 2a ^ —ZtiMxh (m 2 , m 2 L ,m 2 fL ,m 2 iR 



where I = jx or e. Diagram 1(b) provides a contribution given by 

oil , \ 

e 2b ~ -—ZyMih {/,m\ L ,m\ L ,Ml) . (9) 

Diagram 1(c) yields 

e 2c ~ ^iliMih (/, ml t , mj T , Ml) . (10) 
Finally, the contribution of diagram 1(d) is found to be 



e 2d ~ ^LcM 2 f 2 {/,m\ L ,m\ h ,Ml) . (11) 

In these equations, Mi )2 are the U(l) and SU(2) gaugino masses, £ is defined in Eq. (|5]), 
and the function f 2 is defined such that 

ft U J\ - QlQ g( Q ) , frlOg(b) . / 7 7\ / 10 N 

- f 2 {a,b,c,d) = — - -77 + 77 777 777 + {a ^ c,b ^ d). (12) 

(a — b){a — c)[a — a) (b — a)(b — c)(b — a) 



1 There are additional terms which can be written, but these are either L-conserving or subleading in 
the LFV calculation that follows. 



4 



The function f 2 is positive definite and we note several interesting limits in its behav- 
ior. When a = b = c = d, f 2 (a,a,a,a) = l/(6a 2 ); and when a ^> b = c = d, then 
f 2 (a,b,b,b) = \/{2ab). The total contribution to e 2 is the sum of the individual pieces 
listed above. Whether we identify the slepton £ = jl or e depends on the decay process 
we are considering. 

It is clear from the Lagrangian in Eq. (|7]) that the charged lepton masses cannot be 
diagonalized in the same basis as their Higgs couplings. This will allow neutral Higgs 
bosons to mediate LFV processes with rates proportional to e 2 . But in order to proceed, 
we will choose a specific process to study, namely r — > 3/x. Our discussions will generalize 
to other related processes (such as r — > jiee) very easily. 

We will only mention in passing the contributions to e 1; since they do not induce LFV. 
The diagrams which contribute to e\ are (mostly) those of Fig. 1 without the slepton mass 
insertion. The contributions of these diagrams to e\ are found to be [pX 



where 



Cl = —fiM, [2h (Ml mf L , mjj - h (M 2 , /i 2 , mfj + 2j x (M 2 , /i 2 , mj R 
+^M 2 [h (^ 2 , m} L , Ml) + 2h (fi 2 , ml, M. 



- fi(a,b,c) = 



ab\og(a/b) + bc\og(b/c) + ca log(c/a) 



(13) 



(14) 



(a — b)(b — c)(c — a) 

These terms will generate a mass shift for the charged leptons that will appear in our 
final formulae second-order effect. 



Flavor— Violating Tau Decays. We begin by extracting the Tr [il terms in the 
effective Lagrangian of Eq. (|7]). The algebra is similar to that in Ref. [JTTJ] and so we skip 
directly to the result. The relevant LFV interaction has the form: 



C ~ {2G 



,V4 



m T K 32 
cos 2 (3 



cos(/3 - a)h° - sm{(3 - a)H° - iA°] + h.c. (15) 



where 



Kij 



(2 



l + (ei + e 2 {YjY v ) 33 ) tan/? 



(yX 



(16) 



(The Lagrangian for (rRe^)-Higgs can be derived from this by replacing k 32 with k 3 i.) 

Given an LFV (r^ei)-Higgs interaction, the decay r — > 3fi can be generated via 
exchange of h°, H° and A as in Fig. 2. The diagram in Fig. 2 is straightforward to 
calculate. We derive a branching fraction for r — > 3/z of 



Br(r -> 3//) 



-K 



82 



tan 6 /3 



(17) 



768n 3 m\ 

where r T is the r lifetime. To derive this formula, we took the large limit in which 
a -> /3 — 7r/2. 
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Figure 2: The Feynman diagram contributing to r — > 3/i. The shaded interaction vertex is the 



new vertex derived in Eq. 15 



How large can this branching fraction be? Consider the case in which // = Mi 
M 2 = mj = mp, M R = 10 14 GeV and (YjY u ) 32 = 1. Then e 2 -4x 10" 4 and 




Br(r -> ~(lx 10" 
which puts it into the regime that is experimentally accessible at B-factories over the next 



few years. At LHC and SuperKEKB, limits in the region of 10~ 9 should be achievable [14 
allowing a deeper probe into the parameter space. We can also do better if /i 3> M lj2 — 
m|. Then the bino contribution is enhanced by a factor ///Mi; for M\ ~ 100 GeV and 
fi ~ 1 TeV, one can get e 2 ~8x 10~ 4 , resulting in a branching fraction 4 times that 
stated above. However we note that the value of e 2 — 4 x 10 -4 is remarkably stable to 
changes in the SUSY spectrum apart from this large-/i option. 

Discussion. We have demonstrated that LFV in the sleptons can generate large 
LFV in the couplings of leptons to neutral Higgs bosons. It is already well-known that 
sleptonic flavor violation can induce LFV in certain magnetic moment transitions such 
as t — > fij, so it is useful to compare this to the r — > 3/i decay which we just found. We 
first note that the two decays possess very different decoupling behavior so that either 
one could be large while the other is too small to observe. The effective operator for 
r — > 3/i is dimension-6: (l/m 2 A )TfiJIfi The r — ► fi'y operator is formally dimension-5, but 
chiral symmetry requires an m T insertion, so that the operator is actually dimension- 
6: (m T /M| USY )r(T Atl '/iF Atl/ where M SU sy represents the heaviest mass scale to enter the 
slepton-gaugino loops. If sleptons and gauginos are light and A is heavy, then r — > /ry 
would tend to dominate; in the opposite limit and with tan/5 large, r — >• 3/i would 
dominate. Because of this different decoupling behavior, it is impossible to correlate 
the two decays without choosing a specific model. Turning this around, observation of 
one or both of these decays can provide insight into the fundamental SUSY-breaking 
parameters. 

The presence of the r — > fi^y operator can also lead to r — > 3/i if the photon goes 
off-shell. However, for this operator the relation between the two branching fractions is 
roughly model-independent [|J: 

Br(r -> 3/i) ^ q QQ3 ( n0 Higgs-mediation). (19) 
Br(r^/i7) 
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If a ratio much larger than 0.003 is discovered, then this would be clear evidence that some 
new process is generating the r — > 3/x decay, with Higgs-mediation a leading contender. 

Our calculation so far has also relied on Yukawa matrix ansatze which reproduce 
the mass matrix of Eq. (|6|). But what if we had made a different choice for YjY v ? 
Another popular option would be the inverted mass hierarchy in which the (1,2) and 
(1,3) elements of YjY u would be 0(1) and the remainder 0(e) ||. Such a matrix could 
lead to observable r — > e/x/x. The constraints coming from /x — > are already strong; 
for Msusy ^ 100 GeV one finds that (YjY„) 2 i has to be < 10~ 2 |T§. But in the large 



^susy, large tan/3 limit, this bound is weakened and r — > e/x/x could dominate. 

There is another interesting test of the inverted hierarchy models presented by Higgs 
mediation. The decay /x — > 3e can also proceed by neutral Higgs exchange. Though the 
electron Yukawa coupling is tiny, this is offset by the extreme precision of rare /x-decay 
searches. In particular, we find that for (YjY l/ ) 2 i = 1, 

s , u , /tan/3\ 6 /100GeV\ 4 . . 

Br(/x - 3e) ~(5x lO" 14 ) x f j x j . (20) 

Again, this process can generated by the similar /x — > e7 operator taken off-shell, but 
there the ratio is again model-independent: 

Brf/x 3 e ) 

— Y- '- ~ 0.006 (no Higgs-mediation). (21) 

Br(/x — > e7j 

Any deviation from this fixed ratio would, as for the taus, be a strong indication of new 
physics such as that found here. 

Finally, we mention here that one could also calculate the rate for processes involving 
LFV Higgs couplings at both vertices, though we will leave this computation to a future 
work. For example, a second way to generate r — > e/x/x would be to use the r^/x^-Higgs 
coupling that we have been considering in this paper, along with a /x^e^-Higgs coupling. 
We can also generate r — > ee/x with a Tr e^-Higgs coupling at one end and a /I R eL-Higgs 
coupling at the other. But as above, both processes are constrained by non-observation 
of /x — > e7 since they require large (Am|)2i. It is notable that these processes have a 
remarkable tan 8 f3 dependence; however the additional powers are mitigated by additional 
loop suppressions. 
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